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Abstract 

This study investigates the adjustment of large-scale localized buoyancy 
anomalies in mid-latitude regions and the nonlinear evolution of associated 
condensation patterns in both adiabatic and moist-convective environments. 
This investigation is carried out utilizing the two-layer idealized moist- 
convective thermal rotating shallow water (mcTRSW) model. Our investiga- 
tion reveals that the presence of a circular positive potential temperature 
anomaly in the lower layer initiates an anticyclonic high-pressure rotation, 
accompanied by a negative buoyancy anomaly in the upper layer, resulting in 
an anisotropic northeast-southwest tilted circulation of heat flux. The evolu- 
tion of eddy heat fluxes, such as poleward heat flux, energy, and meridional 
elongation of the buoyancy field, heavily depends on the perturbation's 
strength, size, and vertical structure. The heatwave initiates atmospheric insta- 
bility, leading to precipitation systems such as rain bands and asymmetric 
latent heat release due to moist convection in a diabatic environment. This cre- 
ates a comma cloud pattern in the upper troposphere and a comma-shaped 
buoyancy anomaly in the lower layer, accompanied by the emission of inertia 
gravity waves. The southern and eastern sectors of the buoyancy anomaly 
show an upward flux, generating a stronger cross-equatorial flow and inertia- 
gravity waves in a southward and eastward direction. Furthermore, the simu- 
lations reveal a similar asymmetric pattern of total condensed liquid water 
content distribution, accompanied by the intensification of moist convection 
as rain bands. This intensification is more pronounced in barotropic structures 
than in baroclinic configurations with stagnant upper layers. This study high- 
lights the importance of considering moist convection and its effects on atmo- 
spheric and oceanic flows in mid-latitude regions, as well as the role of 


buoyancy anomalies in generating heatwaves and precipitation patterns. 
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1 | INTRODUCTION 


Heatwaves are becoming increasingly frequent and severe 
worldwide, with mid-latitude regions being particularly 
affected. Recent examples include heatwaves observed 
in North America, Europe, and Asia (Diffenbaugh & 
Ashfaq, 2010; Hansen et al., 2012; Perkins et al., 2012). 
The origins of heatwaves can be attributed to various 
physical drivers that fall into two primary categories: 
atmospheric drivers, which encompass processes related 
to moisture and radiation, and surface drivers. Addition- 
ally, the drivers from both the atmosphere and surface can 
interact and produce combined nonlinear effects, as dis- 
cussed in previous studies (for a comprehensive review of 
heatwave drivers, see Horton et al. (2016), Miralles et al. 
(2019). It is important to emphasize that, in the context of 
this study, the term “heatwaves” refers to the intricate pro- 
gression and dynamic evolution of extensive, localized 
positive buoyancy anomalies. These anomalies originate 
from elevated potential temperatures, emerging in the 
absence of simultaneous background zonal fluid motion, 
through “conceptual” examinations, which involve the 
use of an idealized model (also refer to Section 2). Such 
distinct anomalies demonstrate the ability to persist over 
an extended period spanning multiple days. The associa- 
tion between warm extremes and persistent high-pressure 
anomalies, commonly referred to as blocking-like circula- 
tion patterns, has been frequently documented in the 
scientific literature (Cassou et al., 2005; Meehl & 
Tebaldi, 2004; Stefanon et al., 2012; Xoplaki et al., 2003). 
These blocking patterns typically manifest as ridges or 
cut-off highs in the geopotential height contours of the 
upper troposphere. However, the underlying causal mech- 
anisms linking these phenomena remain inadequately 
explored, with limited explicit ridge simulations or con- 
ceptual frameworks available. To address this knowledge 
gap, the present study endeavors to investigate the poten- 
tial influence of surface warm extremes on the occurrence 
of upper tropospheric atmospheric ridges. Our research 
uncovers that upper tropospheric atmospheric ridges can 
manifest concurrently with surface warm localized anom- 
alies. This concurrent development is attributed to the 
heating effects near the Earth's surface, which induce the 
formation of a thermal low-pressure system and a cyclonic 
circulation pattern in the lower atmosphere, countering 
the prevailing upper-level anticyclone. The simultaneous 
occurrence of these phenomena is corroborated by obser- 
vational studies, such as the comprehensive investigation 
conducted by Pfahl and Wernli (2012). It is important to 
note that a universally accepted theory encompassing all 
aspects of heatwaves is currently lacking. The primary 
objective of our study is not to replicate real heatwaves pre- 
cisely but rather to simulate the concept of buoyancy 
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anomalies in a robust manner to capture their fundamental 
characteristics. 

In a “dry’ environment, the establishment of geo- 
strophic balance involves the development of a horizontal 
temperature gradient due to a localized heat disturbance 
(Holton, 2004). The term “dry” signifies an adiabatic pro- 
cess in which no condensation or phase transition of water 
vapor occurs, thereby ensuring the conservation of total 
energy. However, the nonadiabatic adjustment process is 
also influenced by the vertical structure of the atmosphere, 
air mass stability, and the presence of moist convection 
and wind shear. Moist convection has a significant impact 
on the dynamic response of the mid-latitude atmosphere, 
affecting weather and climate prediction and dynamics. To 
accurately model the interactions between moist convec- 
tion and geostrophic adjustment, sophisticated numerical 
modeling techniques are necessary, along with a compre- 
hensive understanding of the fundamental physical pro- 
cesses involved. 

The rotating shallow water (RSW) model is a widely 
used tool for studying atmospheric and oceanic flows in 
mid-latitude regions. It is based on the integration of atmo- 
spheric primitive equations using pseudo-height isobaric 
vertical coordinates, which allows it to capture important 
features of large-scale fluid dynamics such as jet streams 
and Rossby waves (Zeitlin, 2018). However, the classical 
RSW model has some limitations. In particular, it neglects 
the effects of horizontal gradients of potential temperature 
and moist convection, which can be significant in certain 
atmospheric and oceanic flows. These effects arise from the 
interaction between temperature and water vapor content, 
and can lead to the formation of clouds, precipitation, and 
other atmospheric features. To overcome these limitations, 
researchers have developed modified versions of the RSW 
model that take these effects into account. One such model 
is the moist-convective thermal rotating shallow water 
(mcTRSW) model (cf., Rostami et al., 2022). The mcTRSW 
model incorporates the effects of horizontal gradients of 
potential temperature and moist convection, resulting in 
more accurate simulations of atmospheric currents than 
those produced by classical RSW models. Given these 
advantages, the mcTRSW model is a suitable choice for 
investigating the fully nonlinear adjustment of large-scale 
positive buoyancy anomalies in a moist-convective envi- 
ronment. Positive buoyancy anomalies are regions where 
air parcels are warmer and less dense than their surround- 
ings and thus tend to rise. Understanding the dynamics of 
such anomalies is important for predicting the behavior of 
the atmosphere in various weather and climate conditions. 
The mcTRSW model can help shed light on these processes 
and provide insights into the complex interactions between 
temperature, water vapor content, and large-scale atmo- 
spheric dynamics. 
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Several studies have explored the atmospheric 
responses to sea surface temperature (SST) anomalies, 
including moisture, using complex models, as comprehen- 
sively reviewed by (Zhou, 2019). However, within the 
realm of conceptual models, although extensive research 
has been conducted on simplified models such as one- 
dimensional models (e.g., Blumen, 1972; Gill, 1976, 1982) 
and linearized models specifically designed for “dry” 
environments (e.g., Lim & Chang, 1981; Matsuda & 
Takayamal, 1989), the comprehensive examination of 
fully nonlinear adjustments of large-scale positive buoy- 
ancy anomalies using a multi-layer moist-convective ther- 
mal rotating shallow water (TRSW) conceptual model or 
other similar conceptual models are still relatively limited 
in the existing literature. This study aims to fill this 
research gap by investigating the dynamic behavior of 
such adjustments and their potential implications for 
weather and climate prediction. 


2 | MODEL DESCRIPTION 
AND SETUP OF NUMERICAL 
EXPERIMENTS 


The model utilized in this study is supported by a 
theoretical foundation with a historical record for the 
incorporation of moist convection. The model is based 
on the mcRSW model, which has already incorporated 
phase transitions of water vapor and the related latent 
heat release (Bouchut et al., 2009; Lambaerts et al., 
2012; Lambaerts, Lapeyre, & Zeitlin, 2011; Lambaerts, 
Lapeyre, Zeitlin, & Bouchut, 2011). The evolution of the 
vertically averaged humidity and the corresponding 
effects of condensation and latent heat release on the air 
column are defined in the mcRSW model in a simple 
and consistent way, following the seminal ideas of Gill 
(1982). The mcRSW model has been continuously 
improved, and the applied version in this study includes 
additional features such as precipitable water, vaporiza- 
tion, and precipitation, as proposed by Rostami and Zei- 
tlin (2018). This improved model has been used in 
numerous studies to investigate the effects of moist con- 
vection on the dynamics of large-scale Earth and 
planetary jets and vortices (Lahaye & Zeitlin, 2016; 
Lambaerts et al., 2012; Rostami et al., 2017; Rostami & 
Zeitlin, 2017; Rostami & Zeitlin, 2019a, 2019b, 2020, 
2022a). 

The simplified two-layer mcTRSW model is character- 
ized by the set of equations outlined by Rostami et al. 
(2022) as below: 


dyv, 


1 
ap tLe = Vr — iV (in + ha), (1a) 
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d i 
wr + f2x v= —5I2Vb2 — V (nib) — beV (Its +h) 
_ 
Fp, 2 - WC), 


_ 
d;hy + V-(hyv,) =" (-C) +(1-7*) ; 


by Tr 
(1c) 
1 hy — (By /b2)H 
Ota Vea) re) /’) 2 Pal 2) = 
(14) 
1 b, — (H,/h,)B 
Arby +¥1-Vbi =7-(+€) —— a VB (ae) 
1 r 
1 by — (H2/h2)B 
Oita an Pa | C) 2 _ 2) a (1f) 
Aq, + V-(q,¥1) =Es —C, (1g) 
O1Qo + V- (GoV1) =C. (1h) 


The Lagrangian derivative is given by 
d;/dt=0/dt+v;-V, where v; represents the vertically- 
averaged horizontal velocity in layer i, and V denotes the 
two-dimensional gradient operator. The thickness of each 
layer is denoted by h; and numbered from the bottom to 
the top of the domain. The parameters H; and B; repre- 
sent the thickness and buoyancy at rest state, respec- 
tively. The layer-averaged buoyancy variable, bj, is 
defined as 90; /0s, where 6; is the vertically-averaged hori- 
zontal potential temperature in layer i. The potential tem- 
perature (9) in atmospheric dynamics refers to a 
thermodynamic variable that characterizes the tempera- 
ture of a parcel of air in a way that accounts for the 
influence of pressure changes. It is defined as the tem- 
perature that a parcel of dry air would attain if adiabati- 
cally brought to a reference pressure level, typically 
chosen as the surface pressure. It serves as a conserva- 
tive quantity that remains constant for a dry adiabatic 
process and can be expressed as: 0 = T(P,/P)*/°, where 
T is the actual temperature of the parcel, P, is the refer- 
ence pressure level, P is the pressure of the parcel, 7e is 
the specific gas constant for dry air, and C, is the specific 
heat capacity at constant pressure. g is the acceleration 
due to gravity, and @, is the potential temperature at the 
lowermost layer, typically at the surface. Additionally, C 
and E,; denote the condensed liquid water content 
(CLWC) in the lower layer and the bulk formula for sur- 
face evaporation, respectively. The parameters y and y* 
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represent the contribution of condensation, due to latent 
heat release, and Newtonian cooling to the bulk entropy 
of the air column, which can vary from layer to layer. 
The parameter f denotes the Coriolis parameter, while Z 
represents a unit vector in the vertical direction. 

Note that in Equations (1), we assume bj>> q,/h; to 
simplify the equations, where q, represents the specific 
humidity or water vapor content integrated over the air 
column. However, the magnitude of q,/h; relative to 
other terms, such as bj, is not negligible. Therefore, we 
include q,/h; in the model to accurately capture the effect 
of water vapor. Scaling and parameters are the same as 
(Rostami et al., 2022). The corresponding equations of 
precipitable water, downward draft, Newtonian viscosity, 
and vaporization used in the model but not presented in 
Equations (1) follow Rostami and Zeitlin (2018). 

We parameterize CLWC using the Betts—Miller 
method (Betts & Miller, 1986) and surface evaporation, 
Es, using a standard approach for ocean-atmosphere 
exchanges (Katsaros, 2001). The numerical methods used 
in the experiments follow the Dedalus algorithm, which 
utilizes spin-weighted spherical harmonics (Lecoanet 
et al., 2019; Vasil et al., 2019). 

The thermogeostrophic balance relations in an adia- 
batic environment can be expressed as: 


ee 1 
fox v1 = SHOVE —DOV(AY +h), (2a) 


= 1 
foe x vy) = —hy) Vb)? — Sh} Vb; — by Vhy” 
— pv. (2b) 


Here, the superscript (0) denotes the leading-order 
variables and fy represents the Coriolis parameter at a 
reference latitude. 

To clarify the interplay between positive buoyancy 
anomalies or positive gradients of potential temperature 
and negative pressure anomalies, we present the system of 
evolution equations for a two-layer TRSW model in a “dry” 
environment. The model equations can be expressed as 
follows: 


(0; +v,-V)v. +fz x vi =—(Vp), (3a) 
(0; +2: V)v, +fZ x v2 = —(Vp,), (3b) 
Oh +V- (hiv) =0, (3c) 

Athy + V - (hov2) = 0, (3d) 

0b; + ¥,- Vb; =0, (3e) 
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O,by + V2-Vb2 =0. (3f) 


In the above equations, v; = (u;,v;) represents the hor- 
izontal velocity, and (Vp;) = Vp;—hiVb; represents the 
vertically-averaged horizontal gradient of the pressure 
field. The indices i=1,2 correspond to the lower and 
upper layers, respectively. The expressions for p,, D>, hi, 
and hy are defined as follows: p, =(hj)+hi+hz2)bi, 
Do =hyby +(hp+hy)bo, hy =hy+(1/2)h, +h, and 
hy = hy +(1/2)hp. The quantities (Vp,) and (Vp,) can be 
expressed as  (1/2)hiVbi +b, V(hp +h +h.) and 
V(h1b1) + (1/2)h2Vb2 + b2V (hp +h2), respectively. It is 
important to note that the vector fields of (Vp,;) do not 
correspond to the gradient of a scalar. The evolution 
Equations (3) are satisfied within a defined horizontal 
domain D. The first boundary condition is characterized 
by the absence of normal flux, expressed as: 


n-hwv;=0, x aD, (4) 


where 71 represents the outward unit vector. A suitable 
Hamiltonian, which serves as a comprehensive 
measure of the total energy in the system described 
by Equation (3), encompassing both kinetic energy 
and gravitational potential energy, can be defined as 
follows: 


he Pe ae 
= [ex n, (Svi+hubs +hy (54+h.b.) | ‘ (5) 
D 


The functional derivatives of the Hamiltonian H can 
be obtained as follows: 


6H oH 6H ~ 
Bu, tbe om, Sb, tele (6) 


Here, ¢;=p;+(1/2)v;* represents the Bernoulli head 
in each layer. It is important to note that, in general, 
dh; /Oh, does not vanish when calculating 5H/5h;,. The 
total mechanical energy of the system, as defined 
in Equation (5), remains conserved, ensuring the bound- 
edness of the system's free evolution. The term 
d’x(1/2)hjv? represents the kinetic energy associated 
with an elementary fluid column in the ith layer, calcu- 
lated under the Boussinesq approximation. Furthermore, 
the height of the elementary column above the reference 
level z=0 in the ith layer is denoted as hj, and its abso- 
lute potential energy is given by d’xh;h;b;. This conserva- 
tion applies to the sum of kinetic energy and potential 
energy of the active layers, relative to them being filled 
with stable stratified layers with potential temperatures 
@,>06,. For a detailed proof of this conservation law, 
refer to Ripa (1993). 
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This study employs numerical simulations to investi- 
gate the dissipation of a large-scale, localized positive 
buoyancy anomaly in mid-latitudes under both “dry” and 
moist-convective conditions. To isolate the specific 
impact of the initial potential temperature disturbance on 
the atmospheric response and exclude any potential 
influence from background zonal velocity, we omitted 
the inclusion of background flow and topography in our 
conceptual experiments. This is due to the high sensitiv- 
ity of the mid-latitude atmospheric response to lower 
buoyancy anomalies, particularly over the ocean, to the 
background state (Thomson & Vallis, 2018). By adopting 
this approach, we effectively isolate and elucidate the 
intrinsic processes governing the temporal evolution of 
the buoyancy anomaly, facilitating a direct and focused 
investigation into its intricate dynamics and dissipation 
mechanisms with a high degree of clarity and confidence. 
The initial buoyancy disturbance is expressed by a simpli- 
fied circular a-Gaussian equation (cf., Rostami & Zei- 
tlin, 2017), which is constrained not only in latitude but 
also in longitude. The equation is given by: 


ge /, 1 a 
by=eHy/2e—_I| -4-=|G| = -+2 > $= 1,2, i 
fo o 2 206 2 


= H,=0.5. 


(7) 


Here, b; represents the non-dimensional buoyancy 
anomaly for layer i, e« is the maximum amplitude of the 
anomaly, o is a parameter that determines the width of 
the anomaly and is set to 4, and r is the spherical distance 
from the center of the anomaly to the location of interest, 
calculated using the Haversine formula as 
D(Lato,Lono,Lat,Lon), where Laty=40 and Lony =0. 
The function I(x) denotes the Gamma function and 
G(r,a) is given by G(r,a)=[1/I'(a)| [fe 't* 1dt. The 
parameters € and o determine the amplitude and steep- 
ness of the thickness anomaly, respectively, and H; repre- 
sents the thickness of each layer. In the barotropic 
configuration, by = b,, and in the baroclinic case, by =0. 
The equality in phase of upper and lower-level buoyancy 
anomalies defines barotropic disturbances with symmet- 
ric behavior and horizontal influence. Conversely, the 
inequality in phase, with an absent or reduced upper- 
level buoyancy anomaly, characterizes baroclinic distur- 
bances with asymmetry, and the formation of cyclonic 
and anticyclonic circulation patterns. Two different 
amplitudes are considered for e: 0.1 for a weak anomaly 
and 0.3 for a strong anomaly. In a moist-convective envi- 
ronment, the initial distribution of water vapor in the 
lower layer is chosen to be close to saturation, while 
the upper layer is far from saturation. This setup ensures 
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that condensation primarily occurs in the lower layer. 
The model includes entrainment of CLWC from the 
lower to upper layer through upward convective flux, as 
introduced in Rostami and Zeitlin (2022b). 

The zonal and meridional scales of the circular 
a-Gaussian spatial distribution at half-maximum inten- 
sity are approximately 32°, which is of the same order of 
magnitude as the Rossby deformation radius at midlati- 
tudes in a barotropic atmosphere. This suggests that the 
scale of the initial buoyancy disturbance is comparable to 
the characteristic scale of the deformation field in the 
atmosphere, and thus it can potentially induce Rossby 
waves and trigger other dynamical processes in the 
atmosphere. 


3 | RESULTS 


In the early stages of the positive buoyancy anomaly adjust- 
ment process, a nearly isotropic inflow towards the center 
of the anomaly is observed. This is because the anomaly 
initially creates a region of low pressure that attracts sur- 
rounding air masses from all directions. However, as the 
anomaly develops, wind directions begin to deviate from 
the center, primarily due to the Coriolis force, which 
results from the Earth's rotation. In the Northern Hemi- 
sphere, the Coriolis force deflects the direction of air 
masses to the right, while in the Southern Hemisphere, it 
deflects them to the left. As a result, asymmetric conver- 
gent zones develop, where air masses converge toward the 
center of the anomaly from specific directions, leading to 
the development of asymmetry in the CLWC and diver- 
gence fields, as illustrated in Figures 1, 2. These convergent 
zones influence the distribution of moist-convective fronts, 
leading to the formation of precipitable water patterns, 
such as spiral rain bands, around the anomaly. These pat- 
terns result from the interaction between the converging 
air masses and the surrounding environment, including 
the beta effect. Specifically, the beta effect describes how 
the Coriolis force increases with latitude in a rotating sys- 
tem with a varying angular velocity, such as the atmo- 
sphere. Additionally, the asymmetric flow generated by the 
buoyancy anomaly produces perturbations in the large- 
scale circulation that can affect the development and evolu- 
tion of other weather systems in the region, such as 
cyclones and anticyclones. The evolution of both the 
lower- and upper-layer buoyancy anomaly and fronts are 
primarily affected by geostrophic and secondarily by ageos- 
trophic circulations. Geostrophic circulations arise from 
the balance between the Coriolis force and the pressure 
gradient force, while ageostrophic circulations result from 
the imbalance between these forces. The interaction 
between these circulations can result in the formation of 
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FIGURE 1 The evolution of buoyancy perturbations, CLWC (represented in colors) in the lower troposphere, and the corresponding 
velocity field (represented by arrows) during midlatitude adjustment in a moist-convective environment. The first and second rows show the 
results for barotropic and baroclinic configurations, respectively. The third and fourth rows show the evolution of CLWC (represented in 
colors) in the lower layer and the corresponding baroclinic velocity field (represented by arrows) for barotropic and baroclinic 
configurations. The panels in each row display the field evolution at times of 6, 13, 26, 43, and 65 days, arranged from left to right. 

The colorbars are adjusted panel-wise for improved visibility. 
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FIGURE2_ = The left panel shows snapshots of the divergence field, while the right panel displays precipitable water in the lower layer, 
both during the midlatitude adjustment of a barotropic configuration in a moist-convective environment. The corresponding velocity field is 
represented by arrows. These snapshots were taken on the fifth day of simulation. 


frontal zones, which separate air masses of different tem- findings suggest that the heatwave triggers the formation 
peratures and densities, contributing to the development of — of mesoscale convective systems, which can produce heavy 
precipitation systems and other weather phenomena. Our rainfall and corresponding rain bands, especially in the 
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The left panel displays the evolution of inertia-gravity wave activity, while the right panel shows the condensed liquid water 


content (CLWC) in the lower layer. Both panels show the results for the barotropic and baroclinic configurations, where the upper layer is 


stagnant, during midlatitude adjustment caused by a strong positive buoyancy anomaly at the initial time. 


southern and western sectors near the buoyancy anomaly 
(Figure 1). The formation of mesoscale convective systems 
is influenced by various factors such as the strength and 
size of the buoyancy anomaly, the vertical structure of the 
atmosphere, and the presence of other environmental fac- 
tors such as wind shear and moisture content. 

During the initial stages of midlatitude adjustment, there 
is an abrupt emission of inertia-gravity waves, which can be 
observed in inertia-gravity wave activity diagrams (Figure 3). 
In this context, inertia-gravity wave activity is determined by 
integrating the divergence of the velocity field. In a moist- 
convective environment, the southern sectors exhibit a 
heightened occurrence of vigorous inertia-gravity waves in 
comparison to the northern regions. The convective activity 
linked to the buoyancy anomaly engenders more robust and 
coherent regions of condensation that impact the layer 
thicknesses (h;,h2) and buoyancy variables (b;,b2) as 
described by Equation (1), thereby amplifying the wave 
generation in those directions. When a disturbance with 
a large amplitude is introduced in the mid-latitudes, the 
inertia-gravity waves produced during the geostrophic 
adjustment process disperse equatorward and dissipate 
over the tropical region. These waves play an essential 
role in the redistribution of heat and momentum within 
the lower troposphere. The findings of our study indicate 
that mid-latitude heatwaves can enhance the moisture 
transport from the ocean to the atmosphere, leading to 
increased atmospheric moisture content and the potential 
for increased cloud formation and condensation. The 
dominance of southward-propagating inertia-gravity 
waves, along with westward-travelling Rossby waves, is 
also evident in the spatial distribution of the energy den- 
sity anomaly depicted in Figure 8. 


Intense buoyancy anomalies exert a pivotal influence on 
the initiation of cyclonic structures (Figure 1). Empirical 
investigations, notably those conducted within the Atlantic 
basin, furnish compelling evidence for the conspicuous link- 
age between explosive cyclogenesis and regions characterized 
by pronounced SST gradients (Sanders & Gyakum, 1980). 
The genesis of such surface temperature anomalies can be 
attributed to diverse factors, including the total surface net 
radiation, the depletion of soil moisture, or the prevalence 
of drought conditions, thereby substantially curtailing latent 
cooling and amplifying surface temperature anomalies 
(Fischer et al., 2007; Zaitchik et al., 2006). 

Within the purview of our inquiry, we ascertain that 
buoyancy anomalies in the lower troposphere instigate 
the formation of an upper-tropospheric anticyclonic struc- 
ture (Figure 4). This finding robustly substantiates the 
concurrent generation paradigm, whereby cyclonic circu- 
lation patterns manifest at lower atmospheric levels con- 
comitantly with anticyclonic circulation patterns at upper 
atmospheric levels. The manifestation of the Pacific 
Northwest heatwave in June 2021 provided compelling 
evidence of the temporal relationship between tempera- 
ture anomalies and the subsequent intensification of a 
low-pressure system at sea-level pressure. Notably, 1 week 
prior to the heatwave's onset on 25 June, temperature 
anomalies surpassed a threshold of two standard devia- 
tions, thereby initiating the process of low-pressure sys- 
tem intensification (White et al., 2023). The phenomenon 
of enhanced ridging in the (mid-)upper troposphere mani- 
fests as an elongated region typified by relatively elevated 
atmospheric pressure vis-a-vis its environs. This configu- 
ration engenders more stable and enduring weather pat- 
terns, characterized by predominantly clear and arid 
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FIGURE4 The evolution of buoyancy perturbations (represented in color) in the upper troposphere and the corresponding velocity 


field (represented by arrows) associated with midlatitude adjustment in a moist-convective environment. The first and second rows show the 


results for the barotropic and baroclinic configurations, respectively. The third row presents the evolution of precipitable water as cloud 


shields, the formation of a comma pattern, and its entrainment in the upper layer for the baroclinic configuration with a stagnant upper 


layer. The panels in each row display the field's evolution at different time points, including 6, 13, 26, 43, and 65 days, from left to right. 


The colorbars are adjusted panel-wise for improved visibility. 


conditions, thereby facilitating the advection and entrap- 
ment of warm air masses. Notably, during the later stages 
of the simulation, the primary anticyclone nucleus 
migrates poleward and gradually engenders secondary 
vortices (second row of Figure 4). The recurrence of a per- 
sistent pattern of mid- to upper-tropospheric anticyclonic 
circulation and enhanced ridging has been frequently 
documented in previous observational inquiries (cf., Black 
et al., 2004; Cassou et al., 2005; Meehl & Tebaldi, 2004; 
Pfahl & Wernli, 2012; Xoplaki et al., 2003). Therefore, 
while we acknowledge and respect hypotheses such as the 
one put forth by Black et al. (2004), which posits that the 
initial manifestation of extreme temperature anomalies is 
primarily ascribable to the enduring anticyclonic circula- 
tion anomaly, our study unveils the plausibility of an 
additional causal pathway. Specifically, our findings eluci- 
date that anticyclonic structures can arise concomitantly 
due to low-level positive potential temperature anomalies. 
Thus, while not discounting the hypotheses like the one 
postulated by Black et al. (2004), which posits that the ini- 
tial manifestation of extreme temperature anomalies is 


primarily ascribable to the enduring anticyclonic circula- 
tion anomaly, our study unveils the plausibility of an 
additional causal pathway. Thus, the origin of these anti- 
cyclonic structures differs from those created by mid- 
latitude jet streams. It is worth noting that the meander- 
ing of mid-latitude jet streams, such as the polar front jet 
stream, can also contribute to the formation of cyclonic 
and anticyclonic structures in the atmosphere (Okajima 
et al., 2021). When a meander or wave in the jet stream 
becomes stationary or slows down, it can lead to the per- 
sistence of high-pressure systems over a specific region, 
resulting in blocking patterns. 

The evolution of precipitable water in the upper tropo- 
sphere results in the formation of a comma cloud pattern, 
as seen in Figure 4. This pattern is a result of the interaction 
between the cyclonic flow in the lower troposphere and the 
anticyclonic flow in the upper troposphere, creating a dis- 
tinct comma shape. Observational studies have reported 
such comma cloud patterns (Carlson, 1980). 

Comparative analysis of the barotropic and baroclinic 
configurations, when coupled with a stagnant upper 
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layer, reveals that inertia-gravity wave activity is stronger 
in the barotropic configuration. This is attributed to the 
resulting temperature gradient, which drives the forma- 
tion of an anticyclone in the upper troposphere. The anti- 
cyclone persists even after the cyclone has dissipated in 
the lower troposphere. In both configurations, the warm 
core in the lower troposphere strengthens as the cyclone 
intensifies, accompanied by the development of a cold 
core in the upper troposphere. However, the evolution of 
buoyancy in the upper layer differs between the barocli- 
nic and barotropic scenarios. In the baroclinic case, 
buoyancy increases towards the northwest, while in the 
barotropic scenario, buoyancy increases towards the 
south and southwest (Figure 4). These differences in 
the evolution of buoyancy between the two configura- 
tions have significant implications for the subsequent 
development and evolution of midlatitude cyclones. Over 
time, the anticyclone evolves into the upper-level trough 
and ridge (Figure 4). These numerical results align with 
the crudest feature of the formation of a cold core in the 
upper troposphere observed in other studies (e.g., Chang 
et al., 2019). 

As Rossby waves’ westward phase velocity reduces 
with increasing latitude, a disturbance's phase in lower 
latitudes will precede that in higher latitudes. Therefore, 
following the development of cyclonic circulation in the 
lower layer, the northeast-southwest tilt of the buoyancy 
disturbance, as observed during evolution (Figure 1, sec- 
ond panel at Time = 13 days), appears. It is worth 
reminding that low-frequency Rossby waves evolve 
slowly. Subsequently, the isentropic surfaces of both 
weak and strong buoyancy anomalies develop a comma- 
shaped structure and are displaced towards the northwest 
by the beta-gyres. The beta-gyre pattern is significantly 
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affected by various factors such as the anomaly’s strength, 
the presence of moist convection, and the vertical struc- 
ture's barotropic or baroclinic nature (refer to Figures 1, 5). 
Enhanced asymmetry, resulting from beta-gyres, leads to 
instabilities developing in the system's core and the forma- 
tion of smaller rotating vortices known as secondary 
eyewall mesovortices. These mesovortices have been fre- 
quently observed in tropical cyclones, as documented in 
previous studies (e.g., Kossin et al., 2002). The meridional 
displacement of anomalies increases with buoyancy inten- 
sity, and strong northwestward-drifting anomalies eventu- 
ally veer eastward. Eastward veering is more prominent in 
barotropic configurations than in baroclinic ones (refer to 
Figures 1, 5). In an adiabatic environment devoid of moist 
convection, the dissipation of buoyancy anomalies exhibits 
a delayed occurrence compared to their counterparts in a 
moist-convective environment. This temporal discrepancy 
is demonstrated in Figures 1, 5, which illustrate the evolu- 
tion of buoyancy anomalies in both moist-convective and 
dry conditions, respectively. 

During the late stages of strong buoyancy anomaly 
evolution, particularly in barotropic configurations, a 
breaking process occurs, which can result in the trans- 
port of a parallel heat flux towards the north. This com- 
plex process is evident in the Hovmdller diagram, which 
illustrates the anomaly's time evolution as it propagates 
northward (see Figure 6). The breaking of the anomaly 
depends on several factors, such as the strength of the 
anomaly, the presence of moist convection, and the baro- 
clinicity of the vertical structure. In barotropic configura- 
tions, the breaking of the anomaly is more pronounced, 
resulting in a more significant transport of heat towards 
the north. This breaking process is thought to be related 
to the development of strong shear zones near the 
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The evolution of buoyancy perturbations (represented in colors) in the lower layer and the corresponding velocity field 


(represented by arrows) during midlatitude adjustment in “dry” environment. The top and bottom rows show the results for barotropic and 


baroclinic configurations, respectively. The panels in each row display the evolution of the fields at 6, 13, 26, 43, and 65 days after the start of 


the simulation, arranged from left to right. 
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FIGURE6_  Hovmiller diagrams were employed to visualize the spatiotemporal evolution of buoyancy anomalies in the lower layer 


under two distinct atmospheric configurations: barotropic (left panel) and baroclinic (right panel). In the baroclinic setup, characterized by a 


quiescent upper layer, the adjustment process takes place at midlatitudes, prompted by the introduction of a positive buoyancy anomaly in a 


moist-convective environment. 
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FIGURE7_ Anasymmetric distribution of integrated CLWC was observed over 2 months for two different configurations: barotropic 


(left panel) and baroclinic (right panel). Both configurations were initialized with a strong buoyancy anomaly, with the red circle indicating 


the location of the large-scale, localized positive buoyancy anomaly at the initial time. 


periphery of the anomaly, which can trigger the forma- 
tion of small-scale vortices contributing to the heat 
transport. 

Buoyancy forcing at mid-latitudes can trigger the 
formation of convective systems, which can lead to the 
formation of rain bands. Positive buoyancy anomalies in 
the lower troposphere over the ocean can arise from 
the warming of the ocean surface, subsequently influenc- 
ing the lower atmosphere. The underlying causes of 
abrupt increases in SST encompass a range of factors, 
such as changes in the transport of heat by the ocean 
(Oliver et al., 2017; Pearce & Feng, 2013), including inten- 
sified boundary currents and reduced rates of heat loss 
from the ocean to the atmosphere, persistent large-scale 
atmospheric synoptic systems, and coupled air-sea feed- 
back processes such as ENSO events (Bond et al., 2015; 
Lee et al., 2010; Rodrigues et al., 2019). This elevated 


ocean surface temperature leads to the transfer of heat to 
the lower atmosphere through processes such as turbulent 
heat fluxes and sensible heat exchange (Perkins, 2015). 
This heating creates a thermal low-pressure system that 
can draw in moist air from surrounding regions. As the 
moist air rises, it can cool and form convective clouds that 
produce CLWC and precipitable waterfronts. The result- 
ing rain bands can bring relief to the high temperatures 
associated with heatwaves. 

During the adjustment of a positive buoyancy anomaly, 
heatwaves can occur, leading to atmospheric instability 
and the formation of precipitation systems, such as rain 
bands. This instability arises from the release of free energy 
as latent heat during the process of moist convection, 
which is triggered by the presence of positive buoyancy 
anomalies. The resulting heating enhances instability and 
promotes further convection, leading to a positive feedback 
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loop. Figure 7 illustrates the formation of rain bands as the 
heatwave propagates. The distribution of total CLWC 
reveals an area of increased moist convection along the 
southern side of the buoyancy anomaly center. This leads 
to a westward veering of the Rossby wave group because 
the lower meridional modes move more quickly westward 
than the higher meridional modes. The nonuniform heat- 
ing distribution from moist convection initiates pro- 
nounced Rossby wave group development, particularly in 
the barotropic configuration, which exhibits a northeast- 
southwest tilted distribution of total CLWC. This develop- 
ment results not only in the propagation of rain bands to 
the west but also in the formation of mesoscale convective 
systems. The adjustment process persists until the pertur- 
bation has dissipated, and the atmosphere has attained a 
state of enhanced stability. The spatial distribution of the 
integrated energy density anomaly in the lower layer 
(Figure 8) provides compelling evidence that the temporal 
evolution of the buoyancy anomaly in the lower tropo- 
sphere is accompanied by the concurrent westward propa- 
gation of Rossby waves and the generation of inertia 
gravity waves, particularly in the southern sector of the 
buoyancy disturbance. These wave phenomena exhibit a 
more pronounced manifestation within the barotropic con- 
figuration compared to the baroclinic configuration. As 
Rossby waves propagate and interact with the buoyancy 
anomaly, they can enhance or deplete the energy density 
depending on their phase and intensity. For example, 
regions, where Rossby waves reinforce the buoyancy 
anomaly, can exhibit higher energy density, while regions 
where Rossby waves counteract or disrupt the buoyancy 
anomaly may display lower energy density. 
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4 | CONCLUSIONS 


The nonlinear evolution of a large-scale localized positive 
buoyancy anomaly at mid-latitudes and the correspond- 
ing condensation patterns within a complete sphere were 
investigated using an idealized two-layer TRSW model in 
both “dry” and moist-convective environments, encom- 
passing barotropic and baroclinic configurations. Our 
study yielded some insights into the intricate dynamics 
and thermodynamic processes associated with the devel- 
opment and propagation of buoyancy anomalies. The 
results show that intense buoyancy anomalies emerge as 
fundamental drivers in the initiation of cyclonic struc- 
tures characterized by thermal lows in the lower tropo- 
sphere. These structures act as catalysts for the 
subsequent formation of upper-tropospheric anticyclonic 
structures. This highlights the significance of atmo- 
spheric ridges in the middle to upper troposphere as 
potential superior predictors of surface warm extremes, 
surpassing the predictive capacity of surface anticyclones. 
During the early stages of positive buoyancy anomaly 
adjustment, a nearly isotropic inflow towards the center 
of the anomaly is observed, exerting an attractive force 
on surrounding air masses from all directions. This initial 
isotropic flow pattern lays the foundation for subsequent 
asymmetrical flow patterns. As the buoyancy anomaly 
evolves, the Coriolis force begins to influence the wind 
direction, causing deviations from the center. This leads 
to the development of distinct asymmetric convergent 
zones where air masses converge towards the center of 
the anomaly from specific directions. The emergence 
of these asymmetrical flow patterns alters the large-scale 


FIGURE 8 


x10 


An asymmetric spatial distribution of the integrated energy density anomaly in the lower layer within a moist-convective 


environment was observed over a monthly duration for two distinct configurations: barotropic (left panel) and baroclinic (right panel). In 


both configurations, a strong buoyancy anomaly was utilized, and the initial temporal snapshot is denoted by the red circle, signifying the 


location of the large-scale, localized positive buoyancy anomaly. 
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circulation, imparting perturbations that can influence 
the development and evolution of other weather systems 
within the region. 

In addition to their impact on local weather patterns, 
our findings indicate that buoyancy anomalies also con- 
tribute to the formation of mesoscale convective systems 
triggered by heatwaves in a moist-convective environ- 
ment. These convective systems can generate heavy rain- 
fall and rain bands, particularly in the southern and 
western sectors adjacent to the buoyancy anomaly. 
Inertia-gravity waves play a vital role in the redistribution 
of heat and momentum within the lower troposphere, 
with their emission occurring during the initial stages of 
midlatitude adjustment. These waves contribute to the 
complex dynamics of the system, influencing the trans- 
port and distribution of energy within the atmosphere. 
The evolution of precipitable water in the upper tropo- 
sphere is a significant factor in the overall behavior of the 
system, resulting in the formation of a characteristic 
comma cloud pattern. This pattern serves as a visible 
manifestation of the intricate interplay between moisture 
dynamics and buoyancy anomalies. 

Comparing the barotropic and baroclinic configura- 
tions, it is observed that inertia-gravity wave activity is 
more pronounced in the barotropic configuration than 
that of baroclinic when coupled with a stagnant upper 
layer. This highlights the sensitivity of wave behavior to 
the specific atmospheric conditions and configuration. 
Both weak and strong buoyancy anomalies exhibit a 
characteristic comma-shaped structure and are propelled 
northwestward by the beta-gyres. This northwestward 
motion contributes to the transport of heat and energy, 
influencing the regional climate dynamics. During the 
late stages of strong buoyancy anomaly evolution, partic- 
ularly in barotropic configurations, a breaking process 
occurs. This breaking process can lead to the transport of 
a parallel heat flux toward the north, contributing to the 
redistribution of energy within the system. Buoyancy 
forcing at mid-latitudes acts as a triggering mechanism 
for the formation of convective systems, which in turn 
can give rise to the formation of rain bands. These rain 
bands contribute to the regional precipitation patterns 
and play a role in shaping the overall weather conditions. 

Given the conceptual nature of this study, it is imper- 
ative to acknowledge the absence of certain factors, 
namely, the background gradient of potential tempera- 
ture, zonal velocity, radiative transfer, and topography. 
These elements will be the subject of investigation in 
future studies employing Aeolus 2.0, an atmospheric 
model renowned for its intermediate complexity. It is 
noteworthy that the mcTRSW model serves as the 
dynamic core within the Aeolus 2.0 framework. By incor- 
porating these aforementioned factors into subsequent 
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research endeavors, our objective is to attain a more com- 
prehensive and realistic understanding of their intricate 
influences on the dynamics and behavior of the atmo- 
spheric system. 
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